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Ch Abstract 

The normal and tangential Casimir force for the rack gear is cal- 
culated numerically in the case of ideal boundary conditions for the 
electromagnetic field — perfect reflection on the boundaries. The re- 
sulting tangential force appears to be essentially shape dependent. 
l/"~) Relatively small shape variations lead to the essential changes in tan- 

C ■ gential force, whereas normal force remains almost unchanged. 

in 

5 1 Introduction 

m 

During last two decades a large number of studies, both theoretical and ex- 
perimental, deals with the Casimir effect [I]. The growing interest to this 
^ problem is motivated by experimental results [21 El HI E] , that provide rela- 

tively precise confirmation of QFT predictions that lays out of the bounds 
of particle physics. On the other hand, Casimir effect leads to the possibility 
of friction-free nanomechanical devices. Recent studies have aimed to find 
configuration with repulsive Casimir force without dielectric fluid. So the 
Casimir force was calculated for different geometry configurations, especially 
for the configurations with repulsive and tangential forces. Some configura- 
tions are smooth (sphere, cylinder, parabolic cylinder, plate |6j[7J|8]), whereas 
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Figure 1: System parameters and shape of edges: a — period, h — depth of 
profile, s — shift, b — width of the gap. Dotted line W denotes the surface 
in the gap between profiled plates. Case A: flat edge. Case B: smooth edge. 



some configurations are sharp (wedge, cone, knife, needle [HI EJ EU HI]) or 
simply rectangle (flat metallic surfaces with n/2 angles between them). 

It is well-known, that in the classical electrodynamics the precise shape 
of bodies (sharp edges, needles, etc) have a significant influence on the cor- 
responding electromagnetic solutions. So the question arises, whether the 
small changes of shape can change the results for Casimir force? From [TT] 
one can conclude that it is quite possible. For instance, even for a cone with 
finite angle we observe additional singularity, caused by the vertex of the 
cone. It should be stressed, that in [8] the case of sharp knife edge is con- 
sidered as a limiting case of parabolic cylinder, with the smooth dependence 
on the parameter of parabolic curve. Our approach is quite different: we 
preserve the "global" shape properties, and change the geometry only in the 
nearest vicinity of the edges. 

Let's consider the rack gear (see Fig. 1. for geometry details). We have 



two profiled plates, parallel to y axis, period of profiles is a, distance between 
plates is b, and the relative shift of plates is s. This geometry is translational 
invariant along z axis (z axis is orthogonal to picture plane). For simplicity we 
consider the ideal case: material of plates is perfect metal (perfect boundary 
reflection), without any realistic frequency dependence. We also will ignore 
the temperature dependence. So we consider the following question: what 
happens, if we change all n/2 angles of this rectangular geometry by the 
edges? We will investigate two cases: case A - - flat edge, where the n/2 
angle is replaced by two 37r/4 angles, and case B — smooth edge, where 7r/2 
angle is replaced by cylinder of appropriate radius (see Fig. 1). 

From the standard explanation of the tangential Casimir force for this 
geometry one can conclude, that it shouldn't be any dependence on the edge 
shape. Indeed, for fixed gap width b in the case of zero shift (s = 0) the 
average distance between our bodies (plates) is less, then in the case of half- 
period shift (s = a/2). So we should obtain the tangential force, that cannot 
depend on the shape details. However, this conclusion is correct only if we can 
neglect side effects at the edges of profile wells, but for the system considered 
the side effects are significant (see sec. 2). Side effects can essentially depend 
on the shape details, so it is possible to observe shape dependence of the 
tangential Casimir force for our geometry. 

2 Calculation method 

One of the most efficient ways to estimate force between two isolated bodies 
is to calculate energy-momentum tensor for the vacuum state on the surface 
W enclosing one of the bodies 



Vacuum expectation value of the energy-momentum tensor can be expressed 
in terms of euclidian Green function via integral over pure imaginary fre- 
quency (X>(/i) = iji. Originally such approach was suggested in p2] for eval- 
uation van der Waals forces between two bodies. Recently this method has 
been advanced for numerical calculation of the Casimir force [13j El- It is 
interesting to note that the selection of the integration contour u;(/x) = ifi 
is not unique. The employment of different variant of the frequency contour 



were considered in [14]. Here the problem is formulated at real "frequency" 
\i and effective complex dielectric permittivity. 

For the problem in question the surface W separating two plates can be 
chosen as y-z plane (see Fig. 1). After integration on y 
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F n = (0\T u \0)dy, F T = (0|T 12 |0) dy (1) 

Jyo Jyo 

for arbitrary y we obtain the normal and tangential force, corresponding to 
one period a in y-direction and unit length in ^-direction. Then we divide 
these values (for one period) by the period length, and so we get the "density" 
of both forces. 

For further computation euclidian Green function Gij(fi,u,v) should be 
constructed, that is transversal: 

d Ui Gij (/!, u, v) = 0, d v . Gij (/!, u,v) = 0, (2) 

it also satisfy boundary conditions 

TiGij (jl, S, v) = = TjGij (fJL, u, s) , (3) 

where f is tangent vector to the boundary surface at a point s, and, finally, 
it should be solution to equation 

AGij (/i, u, v) - fi 2 Gij (/i, u, v) = Sij (u-v), (4) 

where 

oo 

/- ( kk\ e ik xu-v) 

dk [ Si > ~ li) ~Wf 

— oo 

is transversal delta function. One can easily verify, that 

1 f°° 

(0|A(m)^(^)|0) = / dfiG l:i (fx,u,v). (5) 

7T JO 

Energy-momentum tensor is constructed of the derivatives of the left hand 
side on u a and v@ for v — u, so we can estimate force, provided we solve 
equations for euclidian Green function. Now we should take into account the 
translational invariance on z. We perform Fourier transformation along z 
axis. 



Definitely, one can also construct the series along y axis, because our 
geometry is periodic in y-direction, but straightforward solution of two- 
dimensional equation appears to be more simple and effective. After Fourier 
transformation along z axis we get 

G i:j {fi,u,v) = — / dqexp(iq(u 3 - v 3 ))Gij{q, fi,u 2 ,v 2 ) (6) 

In J 

where u 2 and v 2 are two-dimensional vectors, constructed from the first two 
components of u and v respectively. Finally, for the arbitrary term in the 
expression of energy-momentum tensor for vacuum state we obtain 

(0\d Ua Ai(u)d v gAj(v)\0)\ = - — r / dqd Ua d v / d\i Gij(q,fj,,u 2 , v 2 ) 



(7) 
Here d U3 = iq and d V3 = —iq. For distant bodies renormalization of this 
expression is trivial - - subtraction of Green function G^(q,fi,u 2 ,v 2 ) for 
Minkowski space yields finite expression. 

Therefore, renormalized vacuum expectation values (J5J), (J7|) are expressed 
through the difference G\f l {q 1 [i,u 2 ,v 2 ) = Gij(q, [i,u 2l v 2 ) - G'- j (q, [i,u 2 ,v 2 ), 
which is the solution of the homogeneous Helmholtz equation with the bound- 
ary condition TjGlj n (q, fi,u 2 ,s 2 ) = —TjG®j(q,fi,u 2 ,s 2 ), where s 2 lays on the 
plate surface. One of the most effective approach to solving the problem 
for G^ n (q, fi,u 2 ,v 2 ) is the boundary-element method (BEM) [T5J [TBI H23 • In 
the present work we apply slightly modified BEM: instead of standard spline 
approach for the given boundary element we use polynomial approximation, 
based on surrounding elements. It is similar, but not precisely equal to spline. 
Additionally, we impose more conditions, then the number of elements and 
just minimize discrepancy in resulting overdetermined linear system. Both 
modifications increase precision for given number of elements. 

To estimate errors of this computation scheme, we perform calculations 
for the trivial case h = 0, with well-known explicit analytical result. In fur- 
ther calculations we use the number of point per unit length of the boundary 
that in this trivial case yields relative error about 10~ 4 . We also estimate er- 
ror for our nontrivial boundary by increasing points density and subsequent 
comparison of results. It appears, that nontrivial boundary form increase 
relative error up to 10 -3 . 
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Figure 2: Density of normal force as the function of shift s. Rectangle case 
- solid line, flat edge — dashed line, smooth edge — dotted line. 
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Figure 3: Density of tangential force as the function of shift s. Rectangle 
case - solid line, flat edge — dashed line, smooth edge - dotted line. 

3 Results and discussion 



We use natural system of units ft — c — 1 and choose the geometry param- 
eters a — 2, h — 0.5, b — 1. This choice provide us comparable order of 
magnitude for tangential and normal force density. For the trivial case h = 
the density of normal force will be equal to 7r 2 /240 = 0.0411, and this value 
will be used as a reference point. For plates with profile A the edge lengh 
is \[2L = \/2 ■ 0.08 and in case B radius of the corresponding cylinder is 
R = 0.08. 

From Fig. 2 one can easily find, that value of normal force density f n 
(it should be stressed, that here "density" means value, calculated for one 
period, divided by the length of period) for zero shift s=0 (/ n =0.0124) es- 
sentially differs from the half of our reference value (0.0411/2=0.0205). So 
side effects are significant. At the same time the shape effects for zero shift 



can be neglected, because they are even smaller, then trivial estimation of 
these effects calculated with the proximity force approximation (PFA). Let 
us notice that we do not use the PFA directly because of its low accuracy. 
Instead we multiply the normal force density for plates with edge calculated 
with the PFA by the normalization constant defined as the ratio of forces 

for rectangular geometry JV = jjjfa, where f£ FA = J^ (^ + {b + 2h y ) and 
the value f n was calculated with the use of the method described in previous 
section. For both flat and smooth edge we get approximately / n =0.0121, 
whereas PFA estimation yields N ■ f FFA ~ 0.0115 in case A, where 

f PFA = k 2 f 1 1_ 

Jn 720a 2 \(b + 2h- 2Lf (b + 2hf 

§-2L ( 6 4 ^ 2L ( (b + L) 2 



+3 1 H H 4 - — ' — 1 

o 4 V (b + 2h) 4 ) b 2 {b + 2L) 2 V b{b + 2L) 

and N ■ f FFA ~ 0.0111 in case B, where 
/ 
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cPFA _ 7r 
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(!"2i?) , 2 f dy 



\ b4 { (b + 2^R-y/R 2 -(R-y) 2 



, (l- 2 ^) ,o f d y 



( b + 2h ) i J (b + 2h-2(R- ^/R 2 -(R- y) 2 

For the half-period shift s—1 we obtain almost identical behavior. Side 
effects are significant, whereas shape effects are relatively small. The dif- 
ference between flat edge on the one side, and smooth edge and rectangle 
case on the other side can be explained, if we take into account that average 
distance between bodies is definitely greater for the flat edge (in the case of 
half-period shift). So we can conclude, that for normal force shape effects 
can be neglected, the change in force density appears to be about 2-3%. 

On the contrary, for the density of tangential force we obtain essential 
shape dependence (see Fig. 3). Even absolute value of the difference of 
tangential forces between rectangle case, flat edge case and smooth edge case 
is much greater, then for normal force. Let us remind, that for the flat edge 
we just change one 7r/2 angle to two 37r/4 angles, whereas for the smooth 



edge we have no angles at all. So the dependence on the shape of edge seems 
to be quite reasonable. 

It should be noted, that different dependence on the edge shape for nor- 
mal and tangential force is quite reconcilable with energy reasons. Different 
energy functions for different shapes can have almost identical derivatives on 
the variable x (normal force), and quite different derivatives on the variable 
y (tangential force). Even for normal force we observe different dependence 
on the shift s for different edge shapes (see Fig. 2). 

For both normal and tangential force we obtain the regular dependence 
on the edge size. For example, in the case of flat edge and shift s=0.6 we 
get: 



Edge size 
Density of tan- 
gential force 



_0 

0.00178 



0.04 
0.00153 



0.08 
0.00129 



0.12 
0.00110 



For other values of shift s we observe the same type of dependence. And 
for the case of smooth edge we also obtain this regular dependence. 

4 Conclusion 



Direct numerical computations lead us to the conclusion, that (at least for 
the geometry considered) there is essential dependence of the Casimir force 
on the details of geometry. If we change rectangle structures by the structure 
with edges (flat or smooth), it leads to essential variation of the tangential 
force. It should be noted, that variation of the force value isn't proportional 
to the relative size of edges: maximum edge size we used was 0.14 — about 
15% of typical length in our geometry, whereas tangential force for smooth 
edge appears to be 8 times smaller, then in rectangular case. 

The result observed is definitely pure side effect, but sometimes side ef- 
fects play significant role. It also should be mentioned, that we consider the 
case of perfect metallic surface (perfect mirror). The influence of realistic fre- 
quency dependence on the shape effects is not obvious. Probably, it can mask 
all these effects, but this can be specified only by the direct calculations. As 
for temperature dependence, from [TT] we can conclude, that non-zero tem- 
perature can only amplify the effect observed, because for flat surface and 
sharp edges the temperature dependence appears to be quite different. 
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